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Abstract: Enzymatic biosensors are often used to detect trace levels of some specific 
substance. An alternative methodology is applied for enzymatic assays, in which the 
electrocatalytic kinetic behavior of enzymes is monitored by measuring the faradaic current 
for a variety of substrate and inhibitor concentrations. Here we examine a steady-state and 
pre-steady-state reduction of H2O2 on the horseradish peroxidase electrode. The results 
indicate the substrate-concentration dependence of the steady-state current strictly obeys 
Michaelis-Menten kinetics rules; in other cases there is ambiguity, whereby he 
inhibitor-concentration dependence of the steady-state current has a discontinuity under 
moderate concentration conditions. For pre-steady-state phases, both catalysis and inhibition 
show an abrupt change of the output current. These anomalous phenomena are universal and 
there might be an underlying biochemical or electrochemical rationale. 

Keywords: enzymatic kinetics; inhibition; anomalous behavior; electrochemical biosensor 



1. Introduction 

Enzymatic electrodes are not only used to detect trace levels of specific substances, but also used to 
assay enzyme activities. Peroxidases participate in numerous physiological processes, such as cell wall 
synthesis, plant defence and activation of oxygen. Horseradish peroxidase (HRP) is a representative 
peroxidase; the mechanisms of peroxidase-catalyzed reactions have been widely studied using HRP, 
especially in understanding the biological behavior of the catalyzed oxidation of analytes by H2O2 [1-4]. 
Many H2O2 biosensors have been constructed using HRP due to its ready availability in high purity and 
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low cost [5-7]. The amperometric biosensors with various enzymes are very useful for the determination 
of analytes with a very low concentration [8-15]. Enzymatic sensors based on electrochemical reduction 
of H2O2 [16-18] and oxidation of NADH [19-21] are now available for clinical and environmental 
detection of glucose, pyruvate, hypoxanthine, alcohol and formaldehyde. 

The use of the rotating disk electrode for the study of an electrochemical mechanism was derived by 
Levich [22]. The output currents obtained are elucidated using the Koutecky-Levich formalism and the 
ping-pong kinetic scheme for peroxidase [23-25]. In recent years, a large body of research has focused on 
the mathematical modeling of enzymatic biosensors. Baronas and coworkers [26,27] presented a 
two-dimensional-in-space mathematical model of amperometric biosensors. They developed later a 
model based on non-stationary diffusion equations containing a non-linear term related to the 
Michaelis-Menten kinetics of amperometric enzyme electrodes [28]. Subsequently, they also provided 
another model of the biosensor based on the mixed enzyme kinetics and diffusion limitations in the 
case of substrate inhibition [29], the system of non-linear reaction-diffusion equations [30], and the 
synergistic substrate determination [31]. Andreu and coworkers [32] formulated analytical expressions 
describing the voltammetric response of a reagentless mediated enzyme electrode operated under 
rotating disk conditions. Patre and Sangam [33] built a model based on a diffusion mechanism related to 
Michaelis-Menten kinetics, which can be used in a membrane of the biosensor. Loghambal and 
Rajendran [34] worked out a steady-state non-linear reaction/dif fusion equation. They also developed a 
mathematical model of amperometric and potentiometric biosensor [35]. Recently, Garay and coworkers [36] 
presented a comprehensive numerical treatment of the diffusion and reactions within a sandwich-type 
amperometric biosensor. 

In our work, the electrochemical behavior of HRP enzyme activity is reexamined using H2O2 
reduction on the electrode. This enzymatic biosensor is optimized for the detection of HRP activity, 
including the catalytic and inhibitory kinetic behaviors in steady-state and pre- steady-state phases. 
Actually, only enzymatic activity during steady state follows the Michaelis-Menten formalism, and 
anomalous behaviors could be observed in all other kinetic processes, suggesting that there might be an 
unknown mechanism for the regulation of mass transfer or/and kinetic rate. 

2. Experimental 

2.1. Reagents 

Horseradish peroxidase (EC 1.11.1.7, 300 U/mg) was purchased from Aladdin Reagent Co. Ltd 
(Beijing, China). Bovine serum albumin (BSA) was obtained from NewProbe Chemicals Co. (Beijing, 
China). Hydrogen peroxide (30%, w/v solution) was purchased from Yingdaxigui Reagent Company 
(Tianjin, China) and the concentration of the more diluted hydrogen peroxide solutions prepared from 
30% hydrogen peroxide was determined by titration with potassium permanganate. Phenylhydrazine 
was obtained from Sinopharm Chemical Reagent Co. Ltd (Beijing, China). Phosphate buffer solution 
(PBS, pH 7.0) contained 3.8 M KH 2 P0 4 and 6.2 M Na 2 HP0 4 . All other reagents were of analytical grade. 
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2.2. Electrode Modification 

Prior to modification, the glassy carbon electrode surface was polished with 0.05 mm alumina paste 
on a micro-cloth and subsequently ultrasonically cleaned thoroughly with acetone, NaOH (1:1), HNO3 
(1:1), and doubly distilled water and then dried at room temperature. A mixed solution of 2.0 U/mL 
HRP and 0.5% (w/w) BSA was prepared by dissolving HRP and BSA in 0.2 mL PBS (50 mM, 
pH 7.0). A volume of 3 mL of the mixed solution was then added drop wise on the surface of the 
glassy carbon electrode. The surface area of the electrode is 0.28 cm . The HRP/BSA electrode was 
cross-linked by placing the electrode in a closed vessel contained 25% glutaraldehyde and water vapor 
for 20 min and dried at room temperature for 1 h and stored at 4 °C until use. 

2.3. Apparatus 

Cyclic voltammetry and amperometric experiments were performed with a LabChem 10M 
electrochemical workstation (Tianjin, China) and a conventional three-electrode system. The modified 
glassy carbon electrode served as working electrode, a platinum wire as counter electrode, and a 
saturated calomel electrode (SCE) as reference electrode. A 0.10 M pH 7.0 phosphate buffer solution 
(PBS) was used as supporting electrolyte. A magnetic stirrer (approximately 400 rpm) was employed 
during the amperometric measurements. All experiments were carried out at room temperature. 

2.4. Cyclic Voltammetry Test 

The electrochemical properties of the electrode modified were characterized by cyclic voltammetry. 
Electrochemical experiments were performed in a conventional electrochemical cell containing a three 
electrode arrangement and the potential swept from -0.4 to 0.2 V (vs. SCE) at a sweeping rate of 0. 1 V/s. 
The cyclic voltammetry experiment was performed in 5 mL pH 7.0 PBS. All experiments were 
conducted at 25 °C. 

2.5. Amperometric Assay 

Electrochemical experiments were undertaken in the electrochemical cell containing 25 mL of the 
supporting electrolyte. Amperometric measurements were carried out in a stirred solution (approx. 
400 rpm). The electrodes were immersed into a stirred phosphate buffer solution (pH 6.5) containing 
1.0 mM hydroquinone as the electron mediator and an initial baseline current was recorded. 0.1 mM 
H2O2 solution was added to get an increased steady-state current record. Subsequently, solution of 
increasing phenylhydrazine concentration was added with a micropipette in succession and the current 
decrease was recorded. All experiments were conducted at 25 °C. 

2.6. Correlation Analysis 

Pearson correlation was a measurement of the strength and direction of a linear relationship between 
two groups of data. Statistical significance is determined for a j^-value <0.05 for all tests. Strong 
correlation is determined for an absolute r value >0.7; weak correlation > 0.5. The bi-variate correlation 
analysis was carried out by the R statistical package (version 2.13.0; http://www.r-project.org/) [37]. 
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The linear regression analysis was performed by the online Free Statistics and Forecasting Software 
(version 1.1.23-r6; http://www.wessa.net/slr.wasp). 

3. Results and Discussion 

3.1. The Model 

The following reaction sequence can be proposed for the case of mediated electrocatalytic reduction 
of H 2 0 2 at HRP electrode [38,39]: 

HRP(Fe 3+ ) + H 2 0 2 — Compound(I) + H 2 0 (la) 

Compound(I) + QH, — 2 -*~ Compound(II) + Q (lb) 

Compound(II) + QH 2 ► HRP(Fe 3+ ) + Q ( \ c ) 

Q + 2e — ^ QH 2 ( ld ) 

where HRP is horseradish peroxidase. Fe is the HRP cofactor. Compound (I) (oxidation state +5) and 
Compound(II) (oxidation state +4) are oxidized intermediates. QH 2 and Q represent hydroquinone and 
its oxidized form (benzoquinone), respectively. In these reactions, a two-electron oxidation of the 
ferriheme moiety of HRP enzyme is caused by H2O2. The resulting intermediate Compound (I) 
consists of oxyferryl iron and a porphyrin 71 cation radical. Compound (I) is reduced to Compound (II) 
and then it was further reduced to the original form of HRP (Fe ) by the mediator QH 2 . Q is finally 
reduced back to QH2 by a rapid reaction involving the acceptance of two electrons from the electrode. 

The inhibition of phenylhydrazine to HRP is an anticompetitive inhibition [39]. When the inhibitor 
(phenylhydrazine) is added into the electrolyte solution, it can combine with Compound (I) causing the 
decrease of the current /. The possible mechanism of inhibition could be represented as follows: 

k k 
llRP(Fe 3+ ) • HA — ^ Compound(I) — ^ ... 

+ 

I 



f 

Compound( I) - I 



(2) 



where compound(I) - 1 represents complex of compound(I) with inhibitor. 

After adding substrate or inhibitor dropwise, the current rapidly returns toward equilibrium (steady 
state) with increasing time. For the pre-steady-state part of enzymatic reaction, the current obeys 
exponential decay equation, i.e., 

/ = C + Ax exp(-£ / x) (3) 

where / is the faradaic current at time t. iao is the steady-state current. T is relaxation time, being equal 
to the time constant of an exponential return of a faradaic current to equilibrium after adding substrate 
or inhibitor, which reflects the rate at which excited current return to steady- state current. 
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Suppose that a steady-state current (/oo) is the current when the biosensor system approaches to 
equilibrium as f — ► °°: 

C = Jim/(0 ^ 

where /(f) is the current at time t. In practice, the current at the response time is assumed as /oo. is 
defined as the time when the absolute current slope falls below a given small value (e < 0.0001) [27], i.e., 



tj> = lim< t : 

R ,>0 



1 


di(t) 


m 


dt 



(5) 



A reduction current (/oo) observed on electrode is a combination of the mass transport-limited current 
(/oo,l) and the reaction-rate-limited current (/oo,k) [40], accords with the Koutecky-Levich equation, i.e., 

1 _ 1 1 

For a rotating disk electrode, the mass-transport-limited current (/oo,l) depends on the angular rotation 
velocity (oo) of the electrode and the bulk concentration ([H2O2]) of H2O2, i.e., /oo,l is determined by the 
Levich equation [22,41]: 

C L = 0.62nF^Z) 2/3 v- 1/6 co 1/2 [H 2 0 2 ] (7) 

where n is the number of electrons transferred to the enzyme in one catalytic cycle (n = 2). F is the 
Faraday constant. A is the area of the electrode surface. D is the diffusion coefficient of H2O2. v is the 
kinematic viscosity of water. For a given biosensor that is rotating at a constant angular velocity, /oo,l is 
proportional to [H2O2], Equation (7) reduces to: 

i 00 , L =o 1 [H 2 0 2 ] (8) 

where a x is constant {a x = 0.62nFAD 2/3 v~ l/6 w m ). 

When the transport of H2O2 is high enough to keep its concentration at the electrode surface equal to 
that in the bulk solution, the reaction-rate-limited current /oo,k is governed by the ping-pong kinetic 
scheme, being expressed by [24,39]: 

1 = 1 ( 1 , [I] , k 2 +k 3 1 
C )K nYAT ^[H 2 0 2 ] k Y k 2 k 3 [QH 2 ] ^ 

Where /Is the surface concentration of HRP enzyme. [QH2] is the concentration of hydroquinone. [ I ] is 
the concentration of phenylhydrazine inhibitor. Assuming that = \lnFArk\, b = 1/nFArki, and 
c\ = (fe + fe)/(^Fv4/l2^3[QH2])+l/nFv477 Equation (9) can be also reduced to: 

d =,,2 [ikr* [I]+Ci (io) 

Substituting Equations (8) and (10) into Equation (6), the Koutecky-Levich equation can be reduced to: 

r a3 [ik] +6[I]+c ' 

where #3 = #2+ l/tfi. 
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Combining b[ I ] term with constant c\ in Equation (1 1) in a constant phenylhydrazine concentration, we 
obtain ci = c\ + b[ I ]. The dependence of the steady-state current on the substrate concentration is written as: 



1 



= a. 



1 



- + c 0 



L 3 [H 2 0 2 ] 2 



(12) 



Combining a 3 /[H 2 02] term with constant c\ in Equation (1 1) in a constant H2O2 concentration, we obtain 
C3 = c\ + a3/[H 2 02]. The dependence of the steady-state current on the inhibitor concentration is written as: 



— = b[I] + c 



(13) 



3.2. Electrochemical Biosensor 



Horseradish peroxidase (HRP) and bovine serum albumin (BSA) are immobilized on the surface of a 
glassy carbon (GC) electrode. Hydroquinone/quinone mediator facilitates the electron transfer between 
protein and the electrode. Experiments are performed with H2O2 to ascertain if HRP is indeed 
immobilized on the electrode. 

Figure 1. (a) Cyclic voltammogram of the HRP/BSA biosensor in PBS buffer (pH 7.0): in 
the absence (blue curve) and presence (red curve) of 3.2 mM H2O2. Scan rate = 0.1 V/s. 
(b) Current-time response of the HRP/BSA electrode on 19 successive titrations of H2O2 
(0.1 mM each time) and subsequently 1 1 titrations of phenylhydrazine (0.5 mM each time) 
into the stirring electrolyte solution. Applied potential: -0.2 V. 
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500 1000 1500 2000 2500 
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Figure 1(a) depicts the cyclic voltammetry curves of the proposed biosensor in the presence and 
absence of H2O2. A redox couple is observed in the absence of H2O2. However, an enhancement of the 
peak current is observed in the presence of H2O2. On the other hand, it has been reported that 
phenylhydrazine is able to inhibit or abrogate HRP activity by competing with the normal 
hydroquinone [42,43]. The inhibition effect of phenylhydrazine is measured by using the proposed 
biosensor. 

Upon 19 successive additions of 0.1 mM H2O2 substrate to the electrolyte solution, / increases with 
[H2O2]. Subsequently, upon 11 successive additions of 0.5 mM phenylhydrazine inhibitor to the 
electrolyte solution, / decreases with increasing [I] see Figure 1(b). 

The steep increased (or decreased) in magnitude of the current is observed when H2O2 (or 
phenylhydrazine) is added into the electrolyte solution, and then the current returns to equilibrium 
(steady state following an exponential decay with a certain recovery speed (relaxation time t). The 
values of /oo and x are best computed by non-linear regression of data fitted to Equation (2). These values 
at different substrate concentration are listed in Table SI (Appendix); they at different inhibitor 
concentration in Table S2 (Appendix). 

3.3. Electrochemical Enzyme Assay 

From Figure 2(a), the relation between /«, and [H2O2] is not linear. The calibration curve is valid only 
under low concentration of substrate conditions. More specifically, when [H2O2] < 0.4 mM, the relation 
is very close to a straight line (calibration curve): L, = -7.86(±0.074) x [H 2 0 2 ] -0.025(±0.018) with 
r — -0.999; 7? < 0.0001. This linear correlation is a prerequisite for applying the electrochemical 
biosensors. Thus electrocatalysis reduction to H2O2 with HRP electrode can detect trace levels of H2O2. 
In fact, our result shows the steady-state current follows the Michaelis-Menten kinetic model within a 
broader concentration range. Kinetic constants are determined by fitting the initial rate data to 
Lineweaver-Burk plot Figure 2(b). Thus, the current is regarded as a relative enzyme reaction rate. This 
methodology could monitor in real-time subtle changes in enzyme activity. 

The pre-steady-state time course of the / —> ioo process is fitted to Equation (3) for a single exponential 
decay, in which x is the apparent relaxation time for output current (i.e., the reciprocal of the apparent 
first-order rate constant). It is found, however, that x is not a simple function of the substrate 
concentration. As shown in Figure 3, x is low in both the low- and high-concentration regions; 
a relatively larger x occurs in at a moderate concentration ([H2O2] = 0.8-1.4 mM). This suggests that, 
within this concentration range, the rate at which the equilibrium is approached is slowed down. 

If the kinetic control becomes operative, there is a possible explanation: a "shoulder" appears on the 
curve of £0.5 versus [H 2 0 2 ] in the two-dimensional-in-space biosensor model, as has been proposed in 
reference [27]. Where £0.5 is the half time (to, 5 = 0.693t). Baronas et al. [27] showed that the shoulder 
occurs at [H2O2] near to the Michaelis constant K m . At [H2O2] « K m (x = l/k\K m ) the reaction kinetics 
for [H 2 02]is a zero order throughout the biosensor, whereas for [H2O2] » K m (x= l/£i[H 2 02]) the 
kinetics is a first order throughout. At intermediate values of [EfeC^Jthe kinetics undergoes a transition 
from zero order to first order. This phenomenon is also manifested in inhibition assays, so it is not 
obvious whether agreement between our data and the shoulder of to. 5 is not a coincidence. Inhibition 
kinetics also experiences a transition from zero order to first order, but it is not proof. 
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Figure 2. Effect of the concentrations of H2O2 substrate on the steady-state current, 
(a) Lineweaver-Burk plot, (b) Michaelis-Menten plot. The initial four points are fitted by a 
calibration line (red line) with the slopes of -7.86 and intercept of -0.025 (r = -0.999; 
p< 0.0001). 
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Figure 3. Effect of the concentrations of H2O2 substrate on the relaxation time during 
pre-steady-state phase. 
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3.4. Electrochemical Inhibition Assay 

Phenylhydrazine can inhibit HRP activity [42,43]. In our biosensor, i^ decreased steadily as [I] is 
increased. However, the Lineweaver-Burk plot for the substrate dependence shows that the reaction does 
not follow Michaelis-Menten kinetics, as the curve connecting the experimental points is clearly curved 
Figure 4(a). However, there is a more simple linear relationship between /«, and [ I ]. The dependence of 
the rate of /*> change on [I] is illustrated in Figure 4(b). As shown in the figure, the regression line is 
discontinuous at inhibitor concentration of 3-4 mM. The [I]-/oo relation is broken into two individual 
lines. When [I]<3mM, decrease of i^ along a regression line is given by the equation: 
/oo = -1.92(±0.075)x[I]-9.19(±0.03) with r = -0.999; /?< 0.0001. When [I]>4mM, decline of 
along another line is given by the equation: /oo = -1.41(±0.024)x[I]-8.54(±0.1) with r = -0.999; 
p< 0.0001. 

Figure 4. Effect of the concentrations of phenylhydrazine inhibitor on the steady-state 
current, (a) Lineweaver-Burk plot, (b) The relationship between /*> and [I]. The initial 
seven points are fitted by a straight line with a slope of -1.92 and intercept of -9.19 
(r = -0.999; p < 0.0001). The final four points are fitted by a straight line with a slope of 
-1.41 and intercept of -8.54 (r = -0.999; p< 0.0001). Discontinuity occurs at 
[I] = 3-4 mM. 
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A similar phenomenon occurs in the relation between [I] and x as well. In Figure 5, x increases 
exponentially as the inhibitor concentration increases. However, in the concentration of 3-4 mM, the 
points determined experimentally is lower than x value extrapolated from low concentration and thus the 
curve displays a drop. 

Figure 5. Effect of the concentration of phenylhydrazine inhibitor on the relaxation time 
during pre-steady-state phase. It increased exponentially with increase in concentration. 
A "depression" appears on the curve in the concentration of 3-4 mM. 
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4. Conclusions 

We have examined the electrocatalytic behavior of enzymatic biosensors and possibilities for 
measuring the enzymatic activities. The substrate concentration-dependence of such a steady-state 
current is governed by the Michaelis-Menten equation under the condition of constant diffusion, 
suggesting that the output current reflects the relative rate of the enzyme reaction. For our datasets, we 
find that Michaelis-Menten model does not work well at high inhibitor concentrations and works better 
at low concentrations, and the enzymatic kinetics gives the best results in the absence of inhibitor. 
Deviations from this model cause observed sudden changes in the current occurring at intermediate 
concentrations of inhibitor or substrate. An auto-deceleration of the current occurs when the substrate 
reaches a specific concentration ([H2O2] = 0.7 mM). In contrast, an auto-acceleration occurs when the 
inhibitor reaches another critical concentration ([ I ] = 3 mM). 

Electrochemical enzyme assays provide an extremely sensitive measure but suffer from difficulties in 
interpretation of the data. At relatively higher concentrations of H2O2, the steady-state current follows 
the Michaelis-Menten equation (Figure 2) and the relaxation time is low (Figure 3). These indicate that 
no significant F^C^-induced suicide inactivation of HRP enzyme occurs in our experiment. The 
Levich equation and ping-pong kinetics emphasize that the current should be a monotonically 
increasing (or decreasing) function of the concentration of substrate (or inhibitor) either under diffusion 
limitations or under enzyme kinetics. However, peaks, discontinuities and sudden changes appear in 
Figures 3-5 and their definite mechanism is unknown yet. Also, it is not clear whether the same would 
be true for other enzymes. Electrochemical biosensors do not provide any microscopic information, and 
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even assignment of the relaxation modes (whether it is the mode of a mass transfer or of a chemical 
reaction) is not obvious. 
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Appendix 



Table SI. The steady-state current (/oo) and pre-steady-state relaxation time (t) at different 
substrate concentrations ([H2O2]). 



[H 2 0 2 ] 




X 


Accuracy 




mlYI 


OlA) 


(s) 


X 2 /dof 


r 


0.1 


-0.833 


2.996 


3.01 x 10" 4 


0.990 


0.2 


-1.614 


3.513 


4.52 x 10 -4 


0.992 


0.3 


-2.389 


3.507 


4.08 x 10" 4 


0.992 


0.4 


-3.153 


2.614 


2.20 x 10 -4 


0.993 


0.5 


-3.718 


3.307 


9.32 x 10" 5 


0.995 


0.6 


-4.403 


4.370 


1.87 x 10 -4 


0.985 


0.7 


-4.961 


3.059 


3.17 x 10 -4 


0.979 


0.8 


-5.450 


11.232 


5.94 x 10 -4 


0.962 


0.9 


-5.979 


15.859 


3.71 x 10 -4 


0.981 


1.0 


-6.451 


13.940 


8.40 x 10~ 4 


0.964 
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Table SI. Cont. 



TTT /~V 1 

[H 2 0 2 ] 


^00 


T 


Accuracy 




mM 


(|XA) 


(s) 


X /dof 


r 


1.1 


-7.006 


14.331 


6.21 x 10" 4 


0.977 


1.2 


-7.520 


10.307 


5.10 x 10" 4 


0.980 


1.3 


-8.008 


9.009 


3.47 x 10" 4 


0.984 


1.4 


-8.468 


5.952 


2.27 x 10" 4 


0.988 


1.5 


-8.881 


5.278 


2.24 x 10 -4 


0.986 


1.6 


-9.268 


4.542 


3.17 x 10 -4 


0.973 


1.7 


-9.633 


5.355 


1.49 x 10 -4 


0.988 


1.8 


-9.953 


2.840 


4.20 x 10" 4 


0.927 


1.9 


-10.272 


5.987 


3.87 x 10 -4 


0.954 


Table S2. The steady-state current (z 


a,) and pre-steady-state relaxation time (t) at different 


phenylhydrazine inhibitor concentrations ([ I ]). 






[I] 


^00 


X 


Accuracy 




mM 


(UA) 


(s) 


X 2 /dof 


r 


0.5 


-8.348 


33.777 


1.48 x 10" 3 


0.997 


1.0 


-7.235 


44.277 


6.57 x 10 -4 


0.995 


1.5 


-6.240 


63.354 


5.21 x 10" 4 


0.994 


2.0 


-5.201 


112.63 


3.84 x 10 -4 


0.994 


2.5 


-4.402 


132.71 


4.24 x 10" 4 


0.992 


3.0 


-3.543 


181.32 


2.61 x 10 -4 


0.994 


3.5 


-3.069 


179.08 


3.42 x 10" 4 


0.990 


4.0 


-2.893 


117.37 


2.19 x 10 -4 


0.991 


4.5 


-2.128 


223.28 


1.58 x 10" 4 


0.992 


5.0 


-1.602 


292.14 


1.71 x 10 -4 


0.989 


5.5 


-0.716 


545.95 


1.62 x 10" 4 


0.986 
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